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The prospects of detecting the CP-odd Higgs pseudoscalar (A) in the 
minimal supersymmetric model via its decay into a Z boson and the lighter 
^ ■ CP-even Higgs scalar (h) at the CERN Large Hadron Collider are investi- 



gated. The final state of Z —* l + l~ and h — * bb, may provide a promising 
way to simultaneously detect the A and the h. The compact muon solenoid 
detector performance is adopted for a realistic study of observability. In this 
discovery channel, the masses of the h and the A can be reconstructed. The 
impact of supersymmetric decay modes is considered. 
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I. INTRODUCTION 



The most important experimental goal of the CERN Large Hadron Collider (LHC) is to 
unravel the mechanism of electroweak symmetry breaking. In the Standard Model (SM) of 
electroweak interactions, only one Higgs doublet is required to generate mass for fermions 
as well as gauge bosons. One neutral CP-even Higgs boson (if ) appears after spontaneous 
symmetry breaking. Various extensions of the SM have more complicated Higgs sectors and 
lead to additional physical spin zero fields M . 

The minimal supersymmetric extension of the Standard Model (MSSM) || has two Higgs 
doublets with vacuum expectation values v\ and t>2. After spontaneous symmetry breaking, 
there remain five physical Higgs bosons: a pair of singly charged Higgs bosons H^, two 
neutral CP-even scalars H (heavier) and h (lighter), and a neutral CP-odd pseudoscalar A. 
The Higgs sector is strongly constrained by supersymmetry so that, at the tree level, all Higgs 
boson masses and couplings are determined by just two independent parameters, which are 
commonly chosen to be the mass of the CP-odd pseudoscalar (Ma) and tan/3 = t^/fi- 

At the one loop level, radiative corrections from the t-quark and the b quark Yukawa 
couplings substantially modify the tree level formulae for masses and mixing patterns in 
the Higgs sector. Let us briefly review recent studies on the search for the A at the LHC. 
In most studies, the parameters were selected such that supersymmetric particle (SUSY 
particle) masses were large so that Higgs boson decays to SUSY particles were kinematically 
forbidden. The A — > 77 mode might be observable at the LHC for 120 GeV < Ma < 2mt 
and tan/3 close to one |§- ||. The A does not have tree level couplings with gauge boson 
pairs. The one-loop induced A — > ZZ — > decay mode was found to be a useful 

discovery channel at the LHC only if tan (3 is much less than one and Ma is near the ti 
threshold . If tan (3 is large, the rf decay mode might be a promising discovery channel 
for the A in the MSSM |||8||J. It has been suggested that the A might also be observable 
via its bb decays in a large region of the (Ma, tan/3) plane, provided that sufficient 6-tagging 



capability can be achieved [10,11]. For tan/3 larger than about 10, the A — ► discovery 
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mode at the LHC was found to be observable in a large region of the (Ma, tan (3) parameter 
space pyi^j, even when Higgs decays to SUSY particles are significant [0. This channel 
will provide a good opportunity to reconstruct the masses for neutral Higgs bosons with 
high precision. 

There are, however, regions of parameter space where rates for Higgs boson decays to 
SUSY particles, e.g. the charginos and neutralinos, are substantial or even dominant. While 
these decays reduce the rates for SM signatures, making conventional detection of Higgs 
bosons even more difficult, they also open up a number of new promising modes for Higgs 



detection [|13j. A recent study suggests that if the sleptons are light, it might be possible to 



search for the neutral Higgs bosons via their decays into slepton pairs [|14[ . 

The MSSM Higgs pseudoscalar (A) decays into a Z boson and the lightest Higgs scalar 
(h) with a large branching ratio for tan/3 < 4 and Mz + < Ma < 2mj.fl The final state 
where Z — > l + l~ and h — > t + t~, was found to be a promising channel to simultaneously 



detect the A and the h at the LHC [13] . The h — > bb has a branching fraction about 20 times 
larger than that of h — > rf . It is, therefore, natural to ask whether with a good b-tagging 
efficiency it is possible to search for the A and the h in the decay of A — > Zh with the final 
state of Z — > l + l~ and h — > bb. 

In this paper, we assess the prospects for the simultaneous discovery of the A and the 
h at the LHC via the process A — > Zh — > l + l~bb. The compact muon solenoid detector 
(CMS) performance is adopted for a realistic study of observability.^ A brief summary 
of our preliminary analysis, for the case when Higgs boson decays to SUSY particles are 
kinematically forbidden, has appeared in the CMS Technical Proposal ||. In this complete 
analysis, we have carried out detailed studies for dominant backgrounds, and further, have 



1 The H — > ZA decay is possible for Ma less than about 60 GeV, that is slightly above the 
excluded region of the LEP and within the discovery region for the LEP II. 
2 Recently, similar results for this discovery channel were also found for the ATLAS detector [16]. 
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considered the impact of Higgs decays to SUSY particles. Parton level calculations are 
presented in Section II. Results from more realistic simulations are discussed in Section III. 
We end with a discussion of our results in Section IV. 

II. PARTON LEVEL CALCULATIONS 

The total cross section for the process pp — > A — > Zh — > l + l~bb-\-X is evaluated from the 
cross section a(pp — ■> A + X) multiplied with the branching fractions of A — > Zh, Z — > / + /~ 
and h — > 66. The parton distribution functions of CTEQ2L ]I7[ are chosen to evaluate the 
cross section of pp -> A + X with A 4 = 0.190 GeV and Q 2 = M|. We take M z = 91.187 
GeV, sm 2 6 w = 0.2319, Mvk = M z cos9 w , m b = 4.7 GeV, and m t = 175 GeV. 

Gluon fusion (gg — > A) , via the top quark and the bottom quark triangle loop diagrams, 
is the major source for the Higgs pseudoscalar if tan (3 is less than about 4. The amplitude 
for the process gg — > A is a function of the quark mass squared (m 2 ), M\, and tan/3 |p~8| : 

1. For Ma < mt, the t-loop is almost independent of mt, the ggA coupling can be either 
obtained from the low energy theorem of the axial anomaly ||19| - [PHI or from the exact 
calculation for gg — > A at the limit of mf/M\ ^> 1. 

2. At the threshold of Ma = 2m t , the imaginary part of the amplitude turns on. There- 
fore, the cross section is significantly enhanced when Ma is close to 2m t . 

3. When m 2 q is much less than M\, the amplitude is proportional to m 2 [ln 2 (mg/M\) — 
2z7rln(m 2 ,/M^)]. The t-loop dominates in a large region of tan/3. The cross section is 
almost proportional to cot 2 /? for tan/3 < 10. Only for tan/3 close to m t /mb, can the 
6-loop dominate and the total cross section be enhanced by large tan/3. 

Since the Yukawa coupling of Abb is proportional to tan /3, the production rate of the A 
from bb — > A is enhanced for large values of tan /3. If tan (3 is larger than about 7, the A is 
dominantly produced from 6-quark fusion {bb — > A) |22{| . We have evaluated the cross section 
of the A in pp collisions a(pp — > A + X), with two dominant subprocesses: gg — > A and 
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gg — > Abb. The cross section of gg — > is a good approximation to the 'exact' cross section 
of bb — > A [^2] for less than about 500 GeV. QCD radiative corrections which increase 
the gluon fusion (gg — > A) production cross section by about 50% to 80% for tan/3 ~ 1 
23| , |24|| are not included in our computation for either the signal or the backgrounds. 

We have included complete one loop corrections from both the top and the bottom 
Yukawa interactions to the Higgs masses and couplings using the effective potential |25 - 
|27fl. The formulae presented in Ref. have been employed in our calculations. The 
contributions from the D-terms are usually small |13|J28[| , and therefore, are not included. 
The soft SUSY breaking parameters A t and Ab are taken to be zero. We take rrig ~ mj ~ 
rriq = |//| and consider three sets of parameters: (a) /i = 1000 GeV, such that the Higgs 
boson decays to SUSY particles are kinematically forbidden; (b) |/i| = 500 GeV, such that 
the Higgs boson decays to SUSY particles become significant; and (c) = 300 GeV, such 
that the Higgs boson decays to SUSY particles are large and dominant when tan/3 is less 
than about 10. 

At first, let us take m- g ~ mj ~ nig = fi = 1000 GeV, such that the neutralinos (Z iy i = 
1—4), the charginos (Wj, j = 1,2), the squarks (g), and the sleptons (/), are heavy and the 
Higgs boson decays to SUSY particles are kinematically forbidden. With QCD radiative 
corrections [p9i^0| , the branching fraction of A — > bb is reduced by about a factor of 2. 
The bb mode dominates the Higgs pseudoscalar decays for tan/3 larger than about 4 and 
M A < 2m t . 

For Mz + Mh < Ma < 2m t and tan/3 < 4, the branching fraction of A — > Zh is 
comparable to B(A — > bb) if the SUSY decay modes are kinematically forbidden. The decay 
of A — > Zh can become dominant if tan/3 is less than 2. In this region, the lightest Higgs 
scalar mass is modified mainly by the radiative corrections from the top quark Yukawa 
couplings. The leading radiative corrections to the can be expressed in terms of the 
third generation doublet and singlet squark masses (^q l and mi R ) ||, 

M 2 h = hM A 2 + M z 2 + 5 - ^] (1) 
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where 



f = {(M A 2 - M z 2 ) cos 2/5 + S] 2 + sin 2 2(3{M A 2 + M z 2 ) 2 , (2) 



and 



5 = 2 x ln[(l + + %)]. (3) 

167r 2 MH/ 2 sin 2 /5 LV m 2 A m? ;J v ; 

We have, of course, included complete one loop corrections from both the top and the bottom 
Yukawa interactions in our analysis with the formulae in Ref. [pq| . 

The total cross section for the process pp — > A — > Z/i — > l + l~bb + X is shown in Fig. 
1 as a function of Ma for various values of tan/5. We take m~ g ~ m ; ~ ~ rriq = and 
consider three sets of parameters similar to those chosen in Ref. ||13|| , (a) rriq = fi = 1000 



GeV, such that the Higgs boson decays to SUSY particles are kinematically forbidden; (b) 
rriq = n = 300 GeV, and (c) rriq = —fj, = 300 GeV, such that the Higgs boson decays to 
SUSY particles are large and dominant for tan f3 less than about 10. The signal cross section 
is largest for case (a) where no supersymmetric decays of A are kinematically accessible. For 
the = 300 GeV cases (b) and (c), the branching fraction for SUSY decays to charginos 
and neutralinos are larger for the case of positive /x since W\ and Z 12 are generally lighter 
than when \x < 0. The reduction in the cross section due to the opening of the SUSY modes 
is smaller when tan (5 is large because the b- Yukawa coupling, and hence the decay rate via 
the bb channel, grows with tan j3. We also remark that for fi = 300 GeV and M A ~ 300 GeV, 
and tan/5 < 2, the cross section is further reduced because the decay h — > Z\Z\ become 
accessible, so that the B{h — > bb) is reduced by almost 50%. 

In Fig. 2, we present the total cross section for the process pp — >• A — ■> Zh — > l + l~bb + X 
as a function of tan/5 for various values of Ma- Other parameters are the same as those in 
Fig. 1. Note that the M^ is usually enhanced with a larger tan/5. In Fig. 2(a), the cross 
section for Ma = 200 GeV drops sharply for tan (3 > 4 because M z + M^ becomes larger 
than Ma- In Figs. 2(b) and 2(c), the cross section of M A = 200 GeV is smaller than that 
of Ma = 300 GeV for a large tan (3 because the phase space is suppressed when M A is only 
slightly larger than M z + M^. 
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Fig. 3 shows the contours of total cross section (a) for pp — > A — > Z/i — > l + l~bb + X in 
the (M4, tan /3) plane for a = 1, 10 and 100 fb. The SUSY parameters are the same as those 
in Fig. 1. Also shown as the dashed line is the mass contour for Ma = Mz + Mh- If SUSY 
decays of the A and the h are kinematically forbidden, the cross section can be larger than 
100 fb, for M z + M h < M A < 2m t , and tan (3 less than about 2. 



III. REALISTIC SIMULATIONS 



Recent studies on <fi —>■ bb, <fi = H,h,A |TTJ, have shown that more realistic estimations 



with a detector response could remarkably modify the results of the parton level calculations. 
To make such a comparison some detector model has to be employed. The LHC will have two 
general purpose detectors: the CMS and the ATLAS. Both detectors will provide precise 
muon and electron momentum reconstruction as well as reasonable jet energy resolution. 
We have adopted the CMS detector performance parameters || to estimate the signal and 
backgrounds. Similar results can be expected for the ATLAS detector [fj. 



A. Calculation tools 

The PYTHIA 5.7/JETSET 7.4 generator with the CTEQ2L |7| parton distribution 
functions is used to simulate events at the particle level. To simulate the signal, we use 
PYTHIA to generate the subprocess gg — * A, which is the dominant source of the A for 
low tan {3. But only the kinematics, parton showering and hadronization schemes are used 
directly from PYTHIA. The cross section of pp — > A — > llbb + X is very sensitive to the 
choice of MSSM parameters. To be consistent with the parton level calculations, the results 
from PYTHIA are rescaled to the same cross section and the one-loop corrected M^ is 
incorporated in PYTHIA by hand. 



The output from PYTHIA/ JETSET is processed with the CMS JET program [32] which 
has been developed for the fast simulations of the "realistic" CMS detector response. The 
resolution effects are taken into account with the parametrizations obtained from the detailed 



GEANT |33| simulations. The CMSJET also includes some analysis routines, in particular, 
a set of jet reconstruction algorithms. All events containing e + e~ or /i + /U~ pairs and at least 
two jets are reconstructed and stored for further analysis. 

At the first stage, the minimal kinematical cuts (I) are applied: 

P£ > 10 GeV, 

> 20 GeV, |r^' e | < 2.4, 

E 3 T et > 20 GeV, \rf et \ < 2, 
\Mu — M z \ < 6 GeV. (4) 

In principle, the r] coverage of CMS detector for jets is much better, but for \rf et \ > 2 the 
6-tagging capability rapidly deteriorates. The lepton pair serves as an event trigger. The 
cut on Mu is included to remove the reducible backgrounds from other dilepton sources. 



B. ^-tagging 

The 6-tagging capability of any detector is mainly defined by the tracker system quality. 
Preliminary CMS results indicate that the average 6-tagging efficiency at the level of 40% 
with a mistagging probability of about 2% can be achieved quite easily with the microvertex 
detectors || . The algorithm applied combines common impact parameter and lepton tagging 
techniques. These values look quite reasonable in light of the recent CDF results |33] . There 
is some hope that some more sophisticated 6-tagging algorithms will allow one to achieve 
the 6-tagging efficiency of 60% with a purity of 1% for both the CMS and the ATLAS 
experiments. 

Because of the existing uncertainties we prefer to use here a rather artificial and simplified 
but more general way to take into account 6-tagging. Namely, all jets are matched as a 6-jet 
or a non 6-jet when they are generated. In our analysis, a given jet passing kinematical cuts 
is matched as a 6-jet if there is a 6-quark in a cone AR < 0.3, where 

AR = yjffl* -7] b ) 2 + ((f)i et - (j) b ) 2 . (5) 
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Of course, there is some ambiguity in this procedure, but we have checked that more than 
95% of the 6-jets from the h —>■ bb decays are matched correctly. Our procedure then allows 
us to apply a simplified parameterization of the fe-tagging efficiency to the same event sample 
simply rejecting or accepting a given jet in the event with a certain probability. 

C. Backgrounds and optimal kinematical cuts 

Relevant backgrounds can be separated into two groups. The first one contains irre- 
ducible background processes with a Z and two 6-jets in the final state. The second one 
contains all other potential backgrounds without a Z or two 6-jets in the final state. In our 
analysis, we have applied double 6-tagging to extract the signal since after a preliminary 
analysis it quickly became evident that single 6-tagging is not sufficient. 

Fig. 4 shows the invariant mass distribution for various backgrounds (da j 'dM Ub i) with 
a 6-tagging efficiency of 40% and a purity of 2% after the minimal cuts (I). The dominant 
background is the irreducible process pp — > Zbb + X, provided that the 6-tagging quality 
is high enough. We estimate this process contribution using the gg — > Zbb subprocess 
incorporated in PYTHIA |35| . The total production cross section is about 600 pb at the LHC 
and kinematics are in fairly good agreement with the results of our parton level calculations. 
To estimate the background from pp —>■ Zjj + X, j = jets, when one or both non 6-jets are 
mistagged, we use PYTHIA to generate the pp — > Zj + X processes with p^ > 10 GeV. 
The second jet is generated in the parton shower evolution. To avoid the double counting, 
events with two tagged 6-jets in the final state were eliminated from the Zj sample. The 
background from ti — > llvvbb is at the same level of Z+jets after the minimal cuts and can 
be reduced further by a factor of 3 with the additional cut (II) on the missing Et, 

E™ ss < 40 GeV. (6) 

The potentially dangerous background from ZZ — > llbb is below the level of the dominant 
Zbb by a factor of 10-20, so the peak at Mz in the M b i background distribution is nearly 
invisible. 
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We also estimated the backgrounds from Zee, cebb and bbbb. They are well below the 
uncertainty in our estimations of the Zbb background. The additional isolation cut (III) 

AR(l,jet) > 0.3, (7) 

effectively reduces the backgrounds with the leptons originating from c or b hadrons. The 
c-jets are assumed to be mistagged as the 6-jets with a probability of 9%. 

We have unsuccessfully attempted to improve the significance of the signal using harder 
kinematical cuts. It seems that to optimize significance one has to use as soft cuts on p l T 
and EjT 1 as possible taking into account detector limitations. This conclusion is valid over 
almost the whole mass and tan (3 range except some small region for M A > 330 GeV and 
tan/5 > 1.5, where harder cuts provide some small increase in the significance compared to 
the softest cuts used. 

The last comment concerning Fig. 4 is that all background distributions are peaked 
around M mb ~ 200 GeV. The large background makes the observation for the A very 
difficult in this mass region in spite of the sizeable production cross section of pp — > A + X 
and a significant branching fraction of A — > Zh. 

This discovery channel, A — > Zh — > llbb, is especially attractive because it allows one, 
in principle, to reconstruct the masses for both the A and the h. The optimal experimental 
procedure to extract the signal from the background seems to be an iterative one. At 
the first step, all lljj events passing the cuts I, II, III, and with both 5-tagged jets have 
to be preselected. Then, in the case of known M^, one can directly use the restriction 
\Mjj — Mh\ < AMjj to improve the observability of the A in the Mujj distribution. The 
window size AM^- depends on the detector resolution. For the CMS detector we estimate 
AMjj ~ 12 GeV for this particular kinematical region. It has to be stressed that the 
reconstructed M b ~ b has the systematic downward shift about 5-7 GeV when one uses the 
usual cone jet-finding algorithms. In the case that M h is beyond the reach of LEP II, one 
can use some 'sliding window' procedure. Namely, one should scan the relevant Mjj region 
by small steps, every time placing the middle of the window into the new position, trying to 
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observe some bump in the Mujj distribution. If there is some peak in the M^j distribution, 
one should apply the last step to improve observability for the h (if its mass is unknown) by 
histogramming the Mjj distribution only for events with Mujj in some window around the 
M A . 

In Fig. 5, we illustrate the procedure of extracting the signal for M A = 250, 345 GeV 
and tan/5 =1.5. A signals in the first row are obtained using Mjj mass window shifted 
with respect to the Mh position. The second row presents the Mim distributions provided 
the optimal Mjj position is chosen. And the last row demonstrates the h signal over the 
background provided the events are selected within the proper Mim mass range around the 
bumps found. 

D. Results 

We present the 5cx significance contours for the CMS detector with an integrated lu- 
minosity of L int = 10 and 30 fb _1 in Fig. 6. All the cuts I, II, III, have been applied and 
the sliding window procedure of Fig. 5 has been used to enhance the signal. We show 
results for two sets of the MSSM parameters: (1) m g ~ rrij ~ rriq — ix — 1 TeV and (2) 
m g ~ mi ~ rriq = — ix = 500 GeV. 

We assume here that L int = 30 fb" 1 will be accumulated at the low luminosity regime 
(~ 10 33 cm -2 ) and that pile up effects will be small. A higher luminosity might introduce 
additional complications due to the pile-up effects. However, if the reasonable 6-tagging 
efficiency (40% + 2%) could be achieved at a higher luminosity (about 10 34 cm -2 ), the 
region explorable with L int = 30 fb" 1 and a better 6-tagging efficiency 60% + 1% can also 
be probed with the more modest 6-tagging efficiency (40% + 2%) at L int = 100 fb _1 . 

For m g ~ mj ~ rriq = /i = 1 TeV, the signal of A — > Zh — > l + l~bb can be observable for 
tan (3 < 3 and 180 GeV < Ma < 2m t . For m g ~ ~ rriq = —ix = 500 GeV, this discovery 
channel, A — > Zh — > l + l~bb, can be observable for tan/3 < 1.8 and 180 GeV < Ma < 2m t . 
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IV. CONCLUSIONS 



The CP odd state of the electroweak symmetry breaking sector of the MSSM is very 
difficult to observe at the LHC. We have found that pp — > A — > Zh — > l + l~bb + X may 
provide a promising opportunity to simultaneously search for the A and the h at the LHC. If 
SUSY particles are too heavy to be produced via decays of A, the signal of A — > Zh — > l + l~bb 
can be observable for Mz + Mh < Ma < 2m t and tan (3 < 3. 

The impact of SUSY decay modes on the A — > Zh decay is significant. 

1. For nig ~ rrif ~ niq = —\l = 500 GeV, the signal of A — > Zh — > l + l~bb can be 
observable for tan (3 < 1.8. 

2. For a positive /i with m g ~ mj ~ rrig = fi, the signal of A — > Zh — > l + l~bb might be 
observable only if /i > 600 GeV and tan (3 is close to one. 

3. For nig ~ m ; - ~ rrig = \n\ = 300 GeV, SUSY decay modes reduce the signal to below the 
observable level. In this case, new channels to search for the A become available. The 
most promising discovery mode is the 4/ signal from the decay A — > Z 2 Z 2 — > 4/ + -Zi-Zi 
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FIGURES 

FIG. 1. The total cross section for the process pp — > A — > Z/i — ► l + l~bb + X in fb, as a 
function of Ma, for y^s = 14 TeV, m t = 175 GeV, m g ~ ~ m^, and tan/3 = 1, 3, 10, and 30. 
The soft SUSY breaking parameters A t and Af, are taken to be zero. We take m g ~ mj ~ rriq = \ri\ 
and consider three cases: (a) fx = 1000 GeV, (b) fx = 300 GeV, and (c) -/x = 300 GeV. 

FIG. 2. The total cross section of pp — ► A — > Z/i — > l + l~bb + X in fb, as a function of tan/3, for 
v/s = 14 TeV, mt = 175 GeV, and M4 = 100, 200, 300, and 400 GeV. Three cases are considered: 
(a) rriq = fi = 1000 GeV, (b) rriq = \i = 300 GeV, and (c) rriq = — \i = 300 GeV. Other parameters 
are the same as in FIG. 1. 

FIG. 3. Contours in the (M A , tan.fi) plane, for a(pp -> A -> Zh -> + V) = 1, 10, and 

100 fb. Three cases are considered: (a) rriq = p, = 1000 GeV, (b) rriq = ri = 300 GeV, and (c) 
= — n = 300 GeV. Other parameters are the same as in FIG. 1. 

FIG. 4. Invariant mass distributions (da / dM^) for the dominant SM backgrounds after cuts 
from: pp -> Z + 2 (non-b) jets (Zj); qq -> ZZ -> Zbb (ZZ); gg -► Zbb (Zbb) and 55 it (tt). 
Double 6-tagging with 40% of 6-tagging efficiency and 2% of mistagging is assumed. This figure is 
generated from a simulation with the CMS performance. 

FIG. 5. The illustration for the sliding window procedure to enhance the signal to background 
ratio as described in Sec. III(C) of the text. This figure is generated from a simulation with the 
CMS performance, for rrig ~ mj ~ rriq = ri = 1000 GeV, tan/3 = 1.5 and yfs = 14 TeV. The left 
three figures are for Ma = 250 GeV, and the right ones for M A = 345 GeV. The arrows in the 
second row of figures mark the chosen Mu^ windows. 
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FIG. 6. The 5 a contour in the Ma versus tan/3 plane, for pp — > A — ► Z/i — > l + l~bb + X in 
the MSSM, at the LHC. This figure is generated using a simulation with the CMS performance, 
for nig ~ mj ~ rriq = \fj\ with fx = 1000 GeV (solid) and [i = —500 GeV (dashed). The figure on 
the left corresponds to an integrated luminosity of 10 fb -1 , the one on the right to that of 30 fb -1 . 
Case a) is for 60% of 6-tagging efficiency and 1% mistagging, case b) is for a tagging (mistagging) 
efficiency of 40% (2%). 
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